Stability of metallic stripes in the extended one-band Hubbard model 
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Based on an unrestricted Gutzwiller approximation (GA) we investigate the stripe orientation 
and periodicity in an extended one-band Hubbard model. A negative ratio between next-nearest 
and nearest neighbor hopping t'/t, as appropriate for cuprates, favors partially filled (metallic) 
stripes for both vertical and diagonal configurations. At around optimal doping diagonal stripes, 
site centered (SC) and bond centered (BC) vertical stripes become degenerate suggesting strong 
lateral and orientational fluctuations. We find that within the GA the resulting phase diagram is 
in agreement with experiment whereas it is not in the Hartree-Fock approximation due to a strong 
overestimation of the stripe filling. Results are in agreement with previous calculations within the 
three-band Hubbard model but with the role of SC and BC stripes interchanged. 

PACS numbers: 71.28.+d,71.10.-w,74.72.-h,71.45.1r 



I. INTRODUCTION 

The existence of stripes, i.e. antiphase domain walls 
in an antiferromagnet stabilized by doped holes, is now a 
well established fact in cupra,tei&2££& and nickelat o 7 ' 8 ' 9 
materials. 

Direct evidence for static stripes in cuprates has been 
more clearly established in La2-a;Sr a; Cu04 (LSCO) even- 
tually co-doped with Nd or Eu. However incommen- 
surate inelastic neutron scattering signals (widely inter- 
preted as evidence for dynamical stripe fluctuations) have 
been found in both LSCO and YBa 2 Cu307-,5 (YBCO) 
compounds with a remarkably similar phenomenology. 10 
Even more static charge order has been found in YBCO 
without the need for additions like Ndiii 

In LSCO the inverse stripe spacing grows linearly with 
doping x up to x ~ 1/8 and the orientation rotates from 
vertical (i.e. oriented along the Cu-0 bonds) by 45° to 
diagonal for concentrations lower x « O.OSAi 2 *^ More- 
over, from this linear relation it turns out that in the 
LSCO compounds stripes are characterized by one doped 
hole per two unit cells along the domain wall, that is a 
linear density of added holes v w 0.5 (hereafter called 
'half-filled stripe'). 

Regarding the Bi 2 Sr 2 CaCu20 8 (BiSCCO) family it is 
a very difficult material to perform neutron scattering 
so few experiments exist [see the discussion in Ref . \^ ■ 
Interestingly this study finds a peak much broader than 
the momentum resolution leaving plenty of room for in- 
commensurate effects. The existence of incommensurate 
scattering in BiSCCO is also supported by the experi- 
ments of Ref. While most neutron scattering experi- 
ments focus on the structure of spin excitations, the in- 
homogeneous charge distribution as arising from the for- 
mation of stripes has also been detected by local probes 
like NQRiiSJI and NMR. 18 Due to refinements in the ex- 
perimental technique Haase et ali 1 ^ where even able to 
demonstrate a correlation of charge and density varia- 
tions on short length scales. 



All these experiments suggest that stripes are a com- 
mon phenomena of all cuprates families and therefore 
may be related to superconductivity. 

Concerning the theoretical aspects, stripes were pre- 
dicted by the inhomogcneous Hartree-Fock (HF) approx- 
imation in the three-band Hubbard modeliS, the one- 
band Hubbard22*2I and the t — J model. 22 However these 
pioneering studies predicted insulating stripes with a lin- 
ear density of v = 1 added holes per lattice constant 
along the stripe (instead of v ~ 0.5) which led to an 
early rejection of stripes^ 

Recently we have presented a computation of metal- 
lic mean-field stripes 24 within an unrestricted Gutzwiller 
approximation 2 ^ (GA) applied to the three-band Hub- 
bard model. The behavior of the magnetic incom- 
mensurability e = l/(2d) 3 ' 4 ' 5,6 (d is the distance be- 
tween charged stripes in units of the lattice constant), 
chemical potential 2 ^ 7 ., and transport experiments 2 ^ 2 ^ 
as a function of doping has been explained in a pa- 
rameter free wayi^i Addition of RPA-type fluctuations 
within a recently developed time dependent GA^i has 
explained 3 ^ also the evolution of the optical conductivity 
with doping 3 ^* in a broad frequency range. 

In our previous computations 2 ^*^ we fixed the param- 
eter set of the three band Hubbard model to the values 
derived for the LSCO family within the local density ap- 
proximation (LDA)i^ In this sense we view our previous 
work as the last stage of an ab inito determination of the 
electronic structure. On the contrary we would like to 
explore in the present paper the parameter dependence 
of the results. This kind of investigation is interesting 
not only in order to check the robustness of the results 
but also to obtain hints about the expected trends in the 
behavior among the different cuprates materials. The 
complexity of the three band model, however, makes it 
difficult to perform such study. To achieve our goals we 
adopt here a much simpler Hamiltonian, namely the ex- 
tended one-band Hubbard model. Given the popularity 
of the model it is also interesting to explore, to what ex- 
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tent the physics found in the three band Hubbard model 
can be found in the one band counterpart. In the present 
work we find that for realistic parameters a similar phase 
diagram is found in the one-band and the three-band 
models^ however some important details about the sym- 
metry of stripes differ. 

Regarding the material dependence it has been pro- 
posed recently^ that a one-band model description of 
the various cuprate families essentially differs in the ra- 
tio between next-nearest and nearest neighbor hopping 
t' It. In fact it is found that the transition temperature 
scales with t' '/t ranging from t'/t « —0.15 for the LSCO 
single layer compound up to t'/t w —0.4 for the Tl and 
Hg based materials. We show that this parameter plays 
a key role for the stability and electronic structure of the 
stripes suggesting a possible connection between stripes 
and superconductivity. In particular the higher the t' 
the smaller is the optimum filling v of the stripe (Fig.|HJ). 
Additionally we present a thorough comparison of the 
stability of stripe and polaron textures upon varying t'/t. 

The influence of a next-nearest neighbor hopping term 
on the stripe formation has been previously investigated 
by dynamical mean-field theory (DMFT)2£ and the HF 
approximationi22i22iffi It was found^S^i that increasing 
the ratio \t'/t\ leads to a suppression of static stripes. 
This may indicate a more dynamical character of stripes 
in some systems like Tl and Hg based compounds. 

In case of the HF approximation^ vertical stripe so- 
lutions are only favored for unrealistic small values of 
U/t w 3. ..5 whereas a ratio of U/t m 10 is required 
in order to reproduce the low energy spectrum of the 
three-band model4^ Therefore our investigations here 
are based on the unrestricted Gutzwiller approximation 45 
which provides an excellent variational Ansatz for the 
ground state energy of the Hubbard model. We compare 
our results with the HF approximation an show that HF 
cannot account for the site centered stripe topology as 
relevant in the one-band model. 

In Sec.^we briefly outline our method and present the 
results of our calculation in Sec. lllll Finally we conclude 
the discussion in Sec. IIVI 



II. MODEL AND FORMALISM 

We consider the two-dimensional Hubbard model on a 
square lattice, with hopping restricted to nearest (~ t) 
and next nearest (~ t') neighbors 

H = -* E c l» c ^ C 1 ) 

<ij>,a 

Here c^l destroys (creates) an electron with spin a at site 
i, and n i a = cl a c i a . U is the on-site Hubbard repulsion 
which in the following is fixed to the value U/t = 10 as 



relevant for the cuprates (see e.g. Refs. 44,46). The un- 
restricted Gutzwiller approximation (GA) can be imple- 
mented by either a variational Ansata 4 ^ or the Kotliar- 
Ruckenstein slave-boson scheme^ and results in the fol- 
lowing energy functional 

E GA (p,D) = t J2 (2) 

<ij>,<? 

<<ij>>,a i 

where pij >a — (cj „Cj,o) denotes the single-particle den- 
sity matrix, Di are variational parameters related to the 
double occupancy of the sites and the Gutzwiller hopping 
factors are given by 

z GA _ y/(l - pa + lh ) : -l>,,.„ - Dj) + \jDi{pu-a - Dj) 
\J Pii,(T (1 Pii,<r) 

(3) 

The energy functional Eq. J2Jl has to be minimized with 
respect to the double occupancy parameters D and the 
density matrix p where the latter variation has to be 
constrained to the subspace of Slater determinants. For 
technical aspects of this variational procedure we refer to 
Ref. 

The clusters investigated in the following have dimen- 
sions N = N x x L up to 16 x 16 and we use periodic 
and antiperiodic boundary conditions in order to stabi- 
lize homogeneous (metallic) stripe textures. 



III. RESULTS 

A. Charge and spin profile of stripes 

In Fig. ^ we report the typical charge- and spin profiles 
in a cut perpendicular to the stripe for vertical site (SC) 
and bond (BC) centered stripes as computed with HF 
and GA respectively. All vertical stripe solutions investi- 
gated in this paper are metallic (i.e. homogeneous along 
the stripe) except for integer v and we have suppressed 
1-D instabilities by choosing appropriate boundary con- 
ditions (i.e. periodic or antiperiodic in direction of the 
stripes). 

The total charge shows dips at the central site (sites) 
of the antiphase SC (BC) domain walls (Fig. ^,,d) cor- 
responding to the accumulation of holes. As also found 
in previous studies^^*^ the GA charge modulation is 
much weaker as compared to HF. 

Due to symmetry of the antiphase domain walls the 
spin density is zero along the stripe in case of SC stripes 
(Fig.[IJi,-c). Therefore a large charge reduction is the only 
possibility to reduce the HF double occupancy on these 
sites. On the other hand the double occupancy within 
the GA can be separately tuned by the variational pa- 
rameter D leading to a smaller charge modulation ampli- 
tude with respect to HF. In Fig. ^ we show the smaller 
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FIG. 1: Profile of the charge density (a,d), staggered spin or- 
der parameter (b,e), and double occupancy (c,f) of site cen- 
tered (a-c) and bond centered (d-f) stripes. Results are for 
t'/t = —0.2 and n„ = 4 half-filled stripes on a 16 x 16 lattice 
corresponding to d = 4 and x = 1/8. 



double occupancy at the core stripes sites (like site 2) 
for the GA whereas both HF and GA acquire approxi- 
mately the same value within the AF ordered regions. In 
contrast, the BC stripe solution has finite spin polariza- 
tion on all sites (Fig. so that in this case a large HF 
double occupancy is avoided (Fig. Qf). This has impor- 
tant consequences for the respective stability of BC and 
SC stripes within the two approximations as we will see 
below. 



B. Effect of t' on the band structure 

In order to understand the role of if in the context 
of stripe formation we proceed by investigating the in- 
fluence of this parameter on the electronic structure. 
Within our variational (mean-field) approach the un- 
doped system is characterized by a homogeneous spin- 
density wave state which band structure is character- 
ized by a lower (filled) and an upper (empty) Hubbard 
band, separated by a single particle energy gap A. Due 
to the incorporation of correlation effects already on 
the mean-field level, A GA is significantly smaller than 
within the HF approximation where it is simply given by 
A HF = Um and m denotes the magnetic moment. 

The stripe state at finite doping induces the formation 
of two bands within the Mott-Hubbard gap. These are 
almost flat perpendicular to the stripes but some disper- 
sion develops when the stripes become closer, similar to 
what is found in the 3-band models In FigElwe show 
the dispersion for both bands in the stripe direction (y) 
for various values of t' /t. 

Doped holes occupy states in the lower band (here- 
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FIG. 2: Dispersion of the mid-gap bands along vertical 
d — 4 stripes in a 16 x 16 lattice computed within the GA. 
t'/t — 0, —0.2, —0.5. a) Site centered stripes; b) Bond cen- 
tered stripes 
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FIG. 3: Sketch of the amplitudes $i(fe) along site centered 
stripe sites i for selected k-states as indicated in the band 
structure of Fig. [2^. 1.)= k = (0, 0) of the active band; 
2.)= k = (tt/2,0); 3.)= k = (0,0) of the upper band. 
Empty (full) circles correspond to positive (negative) ampli- 
tudes which value is indicated by the radius. Crosses mark 
$i(fc)=0. 



after called the active band) and a half-filled stripe corre- 
sponds to a half-filled active band. Furtheron the band- 
structure of the BC solution (Fig. Eb) displays a gap 
between both bands which is due to a AF spin-density 
modulation along the BC stripe. For both, SC and BC 
stripes, Fig. |3 reveals a significant broadening of the 
active band upon increasing This broadening has 
a strong influence on the kinetic energy of holes in the 
stripe and therefore also on the density of states at the 
Fermi energy and, as we will show below, on the value of 
the optimum v. 

The broadening can be understood by looking at the 
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single particle wave functions (WF) of selected k-points 
of the active band as shown in Fig. |3| for the case of 
SC stripes. The k = WF is characterized by having 
opposite sign at the core of the stripe than on the ad- 
jacent (AF ordered) legs. As a consequence a positive 
next-nearest neighbor hopping (i.e. negative t', cf. Eq. 

leads to a lowering of this state in first order per- 
turbation theory. By the same argument one can show 
that the zone boundary state of the active band is not 
affected by t' since next-nearest neighbor hopping along 
the stripe always connects to a site i with a node in the 
single particle wave function (the shift in energy in Fig. 
12 is due to the fact that Ep shifts downwards with t 1 ). 
For a sizeable SDW gap A on the adjacent legs of the 
stripe the broadening of the active band can be evalu- 
ated as SB = 36t't 2 /A 2 . Since A HF > A GA as already 
mentioned above, the t' band width renormalization is 
expected to be larger in the GA as compared to HF. 

For completeness we mention that similar considera- 
tions also hold for the upper band in the Mott-Hubbard 
gap. From the structure of the k = (0, 0) wave-function 
(cf. Fig. [21 it turns out that a negative t 1 induces a 
shift of this state to lower energies (in Fig. [3 this shift is 
almost compensated by the change in the chemical poten- 
tial). Correspondingly we find a narrowing of the upper 
band SB = -28t't 2 /A 2 which for large t' leads to the fiat 
structure as observed in Fig. [21 



C. Parametrization of e^; GA vs. HF 

To determine the more stable solution among different 
textures we need to compare the total energies per atom, 
E/N, for the same doping, x = Nh/N. Equivalently we 
can compare the excess energy per hole with respect to 
the undoped antiferromagnet (AF), hereafter the "energy 
per hole" : 

* - E ^>- E " (4) 

where E(Nh) is the total energy of the inhomogeneous 
state (generally a stripe state) with Nh holes counted 
from half-filling and Eaf is the AF solution on the 
same cluster. The minimum of en at low doping (non- 
interacting stripes) has the following meaning: it deter- 
mines the optimum filling that stripes should have in or- 
der to accommodate a given number of holes in a variable 
number of stripes. 

In Fig. 01 we report the energies per hole for different 
configurations in the HF and in the GA approximations 
as a function of the stripe filling v — Nh/(n s L). Here L 
is the linear dimension in the y direction of the cluster 
which coincides with the stripe length, n s denotes the 
number of stripes and Nh is the total number of holes 
(counted from half-filling) in the system. More detailed 
results for vertical SC and diagonal stripes in the GA are 
reported in Fig. 



In order to analyze these results we expand the excess 
energy per unit length as a Taylor series in v. 

m h )-E A , =A+Bu+Cu2 (5) 

The energy per hole in a stripe thus follows the relation 

e h = A/v + B + Cv (6) 

which provides an excellent fit to the data points as 
shown in Figs. ^>,|S1 The minimum of is given by 

^A/C 

and the second derivative of the energy per hole at the 
minimum reads as 

e'> mm ) = ICjcjA. (7) 

Note, however, that the single particle spectrum has a 
gap for v = 1. Therefore the energy has a cusp-like 
singularity at this filling which causes the v = 1 state to 
be the minimum for a wide range of parameters. We can 
take this into account by modifying the above formula 
to: 

v min = min(^VC,l) (8) 

In the dilute limit, corresponding to large stripe separa- 
tion, the coefficients can be interpreted as follows: A is 
the energy cost per unit length to create an antiphase do- 
main wall in the stoichiometric (otherwise AF ordered) 
system. The parameter B can be defined as the chem- 
ical potential to add one hole into the "empty" stripe. 
Finally the parameter C — Ck + Cj is related to the 
kinetic (Ck) and interaction (Cj) energy of the system. 

In order to separate the two contributions Ck, Cj to 
the parameter C we have adopted the following proce- 
dure. For a small value of v we have optimized the 
variational parameters of the GA. Then the energy as 
a function of v has been computed without changing the 
variational parameters, i.e. by solely filling the frozen 
bands. The resulting quadratic contribution to the en- 
ergy gives the C'k parameter and since C is known we 
also obtain C\. Within mean-field the latter is due to 
readjustments of the charges and double occupancies in 
response to the added holes. 

For the parameters under consideration we find that Ci 
is approximately 2 — 3 times larger than Ck ■ However, it 
is important to note that Cj is almost independent from 
the next-nearest neighbor hopping, whereas from Fig. [21 
it is obvious that the kinetic energy of the holes moving 
along the stripe will have a significant dependence on t' . 
To see more explicitly how the active bandwith reflects 
on the kinetic energy one can approximate the active 
band density of states per spin by a rectangular shape 
of height 1/W(t'). Here W(t') represents an effective 
bandwith which is roughly given by the width of the lower 
(active) band in Fig. [21 One finds that the kinetic energy 
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FIG. 4: (a) HF energy per hole eu for vertical stripes as a 
function of the stripe filling v and for various values of t'/t. 
Results in (a) are for 2 d — 5 stripes on a 10 x 10 lattice, 
(b) GA energy per hole for 4 site centered (SC) and bond 
centered (BC) stripes on a 12 x 12 lattice (symbols). The 
lines correspond to fits obtained with Eq. @ 



is quadratic in v and Ck ~ W(t'). Below we discuss in 
detail how t' determines the stability of stripes via this 
mechanism. 

To estimate A we need an estimate of the energy to 
create an empty domain wall. Since this is related to an 
excitation of the undoped system we can use the Heisen- 
berg model with a magnetic interaction J = 4t 2 /U in 
the limit of U >> t. The estimate is easily done in the 
case of BC stripes which have a core with ferromagnetic 
alignment along the domain wall. In the Neel limit one 
obtains: A = J/2. 

The case of SC stripes is more subtle. Indeed for 
SC stripes one has a core without permanent moments 
which, contrary to the BC case, does not have a classical 
analog within a Heisenberg description of the undoped 
state. A non classical state will consist of a spin liquid 
along the core which in any case will give A ~ J. Within 
the HF approximation of the Hubbard model the corre- 
sponding solution has a non-magnetic core with maxi- 
mum double occupancy (1/4). This is at odds with the 
Heisenberg picture of a spin liquid core which, by con- 
struction, implies small double occupancy for the original 
fcrmions. Therefore the HF approximation introduces an 
error in A which in the large U limit is of the order of 
[7/4. This is due to the fact that moment formation is 
the only mechanism within HF to reduce double occu- 
pancy. On the other hand the GA allows an unpolarized 
state with small double occupancy (D ~ t 2 /U 2 ) thanks 
to the additional variational parameters Di as shown in 
Fig. nji,f. Due to this mechanism we obtain A within 
the GA of the correct order A ~ UD ~ J. We see in 
Fig- El that the curves for SC and BC stripes are similar 
in the GA whereas in HF metallic SC stripes are energeti- 
cally unfavorable with respect to BC (and diagonal) ones. 



Clearly the latter faultiness of the HF approach is due to 
the overestimation of the parameter A mentioned above 
[c.f. Eq. [8] which leads to an unphysical large value of 
Vrain ■ For this reason we conclude that HF theory is not 
an adequate approach for the description of SC stripes. 

Within the HF approximation it is known^ that for 
U/t > 3.6 vertical stripes are less stable than diagonal 
ones. Moreover spin polarons become the ground state 
when U/t > 8. Fig. demonstrates that similar results 
hold when one includes next-nearest neighbor hopping. 
As we will show below the situation is quite different in 
the GA approximation which also allows a comparison 
between SC and BC stripes on an equal footing. 

It is seen from Fig. 0Jd that within the Gutzwiller ap- 
proach SC stripes become degenerate with BC for large 
and small filling factors v. Especially it turns out that 
BC empty stripes are slightly lower in energy than SC 
empty stripes. However, SC stripes tend to have larger 
bandwith than BC ones (cf. Fig. [2J which leads to a 
larger C, larger curvatures of vs. v [Eq . J7J)] and 
hence smaller values of v. The difference in curvature 
allows for the two intersections between the e% curves for 
BC and SC solutions in Fig. 0Jd. We will show below 
that the crossing at large v corresponds to optimal dop- 
ing and that in the underdoped case SC stripes are more 
stable than BC ones. The stability of SC solutions with 
respect to BC ones at low doping as well as the increased 
degeneracy for larger concentration is in agreement with 
DMFT calculations in the one-band models Therefore 
it seems to be a robust feature of the model and not 
to depend on the approximation. On the other hand it 
is opposite to what is observed in the 3-band model^ 
where BC stripes constitute the most stable low doping 
structure. 



D. Effect of t' on the stability of stripes 

We now proceed to discuss in detail the role of t' on 
e/j. Due to the dependence of Ck on t' and the con- 
stancy of the parameter A we observe in Fig. [S] a shift 
of optimal stripe filling to smaller values and a steeper 
minimum upon increasing Hence, the larger is \t'\ 
the more stable are partially filled stripes with respect 
to their filled (insulating) counterpart in agreement with 
the estimate Ck ~ W ~ \t'\ and Eq. JSJ. One can un- 
derstand the reason for this behavior by comparing the 
two cases v — 1 (insulating) and v — 1/2 (metallic), re- 
membering that the former structure has twice as much 
holes per stripe. For the sake of simplicity let us ne- 
glect for the moment the interaction part Cj. Then, if 
the cost to create a stripe dominates (A > C or roughly 
speaking J > W), v = 1 stripes are more convenient be- 
cause this minimizes the number of stripes that have to 
be created to accommodate a given number of holes. On 
the other hand for v = 1 the Pauli exclusion principle 
forces to fill the higher kinetic energy states of the active 
band, which for broad bands becomes energetically unfa- 
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vorable. As a consequence it is more convenient for sys- 
tems with broad bands (W > J) to pay the energy J in 
order to create more stripes, thus reducing the filling per 
stripe and avoiding the occupancy of high kinetic energy 
states of the active band. Clearly the relevant parame- 
ter which controls this behavior is J/\t'\ ~ J/W ~ A/C 
with small ratio A/C favoring partially filled stripes in 
agreement with Eq. JSJ (see also Fig. below). We note 
that this is in agreement with our findings in the three- 
band Hubbard model where the role of t' was played by 
the oxygen-oxygen hopping. 24 

The energy of (vertical) metallic stripes has to be com- 
pared with other saddle-points of the GA energy func- 
tional. It should be noted that especially at large dop- 
ing (but depending also on t'/t) there exists a huge va- 
riety of inhomogeneous charge and spin textures which 
are rather close in energy and compete with the vertical 
stripe solution. Among this variety we have selected di- 
agonal stripes and spin polarons which energy per hole 
is reported in Fig. [SJd. This choice is motivated by the 
idea that with increasing kinetic energy the stripes should 
undergoe some kind of melting process which reflects in 
substantial lateral and orientiational fluctuations and an 
eventual decay into individual segments. However, in our 
static approach we can only access some 'snapshot' con- 
figurations of this melting process and demonstrate the 
increased degeneracy in energy of these structures with 
doping. 

Two major points result from Fig. |SJ i) For fixed 
stripe periodicity d vertical solutions are in general lower 
in energy than diagonal onesi-^ However, for larger stripe 
separation (d = 5, 6) and small ratio \t'/t\ < 0.1 one ob- 
serves from Fig. [5] that diagonal and vertical structures 
are rather close in energy, ii) For moderate values of 
\t' /t\ ^ 0.4 spin polarons are always higher in energy 
than (vertical) stripe structures but the polaron solution 
becomes more stable for larger \t' /t\ depending on the 
stripe separation. At low doping where the stripes are 
almost noninteracting (d > 5) the critical value for (ver- 
tical) stripe stability is t'/t « —0.4. In this region of 
parameters the ground state will have contributions from 
both stripes and dissociated holes. We can speculate that 
in this case, going beyond mean-field, stripes will have a 
more dynamical character. In this sense our finding is in 
agreement with previous investigations^^ which found 
a suppression of static stripe formation upon increasing 
\t'/t\. Note that for finite doping the polaronic energy 
naturally depends on the special arrangement of the self- 
trapped charge carriers on the lattice. However, for the 
selected large value of U/t = 10 the spin-polarons are 
quite localized and therefore their energy is almost inde- 
pendent of doping up to hole concentrations x « 0.2. 

Fig. El displays the optimal filling v opt (corresponding 
to the curve minima in Fig. |3J versus the stripe sepa- 
ration for different next-nearest neighbor hoppings. For 
stripe distances d > 4 the vertical stripes are almost non- 
interacting and therefore the filling fraction versus t'/t is 
the same as can be seen from Fig. [H^,. For smaller dis- 
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FIG. 5: GA excess energies per hole for vertical SC (a) and 
diagonal (b) stripes as a function of the stripe filling v and 
for various values of t'/t. The stripe separation (in lattice 
constants) is d = 3,4,5,6 (from top to bottom). We also 
report the energy of the polaron solutions in the right panel. 
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FIG. 6: Optimal stripe filling as a function of periodicity for 
(a) vertical and (b) diagonal stripes. 



tances stripes start to repel, resulting in a shift of 
to higher values (cf. Fig. |2K) - However, the increased 
inter-stripe interaction leads also to a larger dispersion 
perpendicular to the stripes which in the curves of 
Fig. [SI reflects as an increase in curvature (i.e. in the 
'kinetic' constant C) upon going from d = 4 to d = 3 
stripes. As a result one therefore observes a shift of v opt 
to lower values for the d = 3 solution. 

The situation for diagonal stripes is qualitatively sim- 
ilar, i.e. upon increasing \t' /t\ the minimum of shifts 
to lower values of v. Since the kinetic energy in diag- 
onal direction (and thus the parameter Cdiag) is natu- 
rally smaller than along the bonds we find the optimal 
filling at larger values as compared to vertical stripes. 
Moreover, diagonal stripes are more extended so that a 
sizeable inter-stripe interaction is present even for peri- 
odicities d = 5, 6. In Fig. |SJd this reflects as a continuous 
decrease of the eh curves upon increasing d. For the same 
reason v opt continuously increases with d for fixed t'/t as 
can be seen from Fig. ^p. 

Comparing Fig. [SJl with the experimental situation in 
the LSCO high-T c compound it thus turns out that the 
phenomena of half-filled vertical stripes for d > 4 requires 
for the ratio between next-nearest and nearest neighbor 
hopping t'/t w —0.2. For this value Fig. [7| reports the 
corresponding eh curves but now plotted as a function of 
the hole doping. Since the charge modulation of verti- 
cal stripes has a width of about 4-5 lattice constants the 
intcr-stripc interaction is vanishingly small for d > 5 and 
thus the energy depends only weakly on distance. More- 
over, from Fig. [3] we have seen that in this regime the 
ground state solution has always filling close to v = 1/2 
and thus the incommensurability 



x 

2^' 



(9) 



behaves as e m x in agreement with experiment For 
d < 4 (x > 1/8) stripes overlap (see Fig.^l which leads to 
a vertical shift of the curves as shown in Fig. EJi- Hence 
it is unfavorable for the system to increase the density 
of stripes from d = 4 to d = 3, the incommensurability 
gets locked at e = 1/8 (d = 4) and the filling v of the 
stripes starts to increase. This feature favors the degen- 
eracy between bond and site centered vertical stripes as 
can be deduced from Fig. 2fc suggesting the importance 
of lateral stripe fluctuations in this doping range. Disre- 
garding for the moment the polaron solution we observe 
a further transition toward d = 3 stripes at x ~ 0.2. A 
similar result has been found also in the 3-band modeli— ■ 
Diagonal stripes are more extended so that the corre- 
sponding eh curves start to shift for d < 5 already. In 
the low doping regime they are less stable than vertical 
stripes, however, the energy difference becomes rather 
small for x > 0.2. Moreover, both vertical and diago- 
nal stripe solutions approach the energy of the polaron 
solution at this particular doping. This means that for 
concentrations x > 0.2 the energy is almost independent 
of the charge topology. The GA landscape is not dom- 
inated by a well defined saddle point and an associated 
Slater determinant but a large number of energy min- 
ima with very different charge topologies compete. In 
this regime the ground state will thus have contributions 
from all Slater determinants corresponding to the dif- 
ferent charge configurations with equivalent energy. We 
expect a state which has lost the orientational order in- 
formation^ 

In the low doping regime neutron scattering experi- 
ments on LSCO observe a change of the incommensu- 
rability vector below x w 0.05 from the vertical to the 
diagonal direction^ii 2 * 1 ^ From the results as shown in 
Fig. [7|we cannot infer a crossover towards diagonal so- 
lutions at low doping. However, the closeness in energy 
of vertical and diagonal structures suggests that the or- 
thorombic distortion, which continuously increases upon 
underdoping in LSCO may lead to a stabilization of the 
diagonal stripe phase at small hole concentrations (see 
Normand et alS for a detailed discussion of the relation 
between stripe orientation and lattice distortion) . More- 
over, we believe that other effects, not included in our 
model, are important in this regard. Indeed the contri- 
bution of long-range Coulomb interactions is expected to 
become important at low doping. A point charge esti- 
mate shows that this tends to favor diagonal configura- 
tions over vertical ones. 

To summarize the results for t'/t — —0.2: We ob- 
tain half-filled vertical stripes and e = x up to hole 
concentrations x = 1/8 and increasing v and e = 1/8 
for 1/8 < x < 0.2 in agreement with experiment 3,4 '-' 6 
and with our previous 3-band computations^ The be- 
havior for the filling implies that the chemical potential 
fi will be approximately constant for x < 1/8 and de- 
crease with doping for x > 1/8 which is also supported 
by experiment. 26,27 However, it should be noted that the 
theoretically obtained shift in \x exceeds the experimen- 
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FIG. 7: Energies per hole for vertical (a) and diagonal (b) 
stripes as a function of doping, t'/t = —0.2 as appropriate 
for LSCO. The inset displays the incommensurability versus 
doping as deduced from (a). The diagonal pattern indicates 
the degeneracy with diagonal stripes. 

tal one by a factor of 2... 3 which may be attributed to 
finite size effects. A similar discrepancy is also observed 
in the 3-band model and we refer to Ref. |2J for a further 
discussion of this point. 

Finally, in Fig. IHlwe show the phase diagram calculated 
for t'/t — —0.4 which is considered to be the appropri- 
ate ratio for the Tl 2 Ba2Ca2Cu30iO and HgBa 2 CaCu206 
compounds. 36 In this case low doping vertical stripes 
have a filling fraction of v « 0.4 so that the incom- 
mensurability e ~ 1.25x is expected to grow faster than 
for the LSCO compounds. However, already the d = 4 
structures are unstable with respect to the spin polaron 
solution and thus the incommensurability starts to satu- 
rate for x > 0.08 corresponding to the d — 5 structure. 
Moreover since the polaron solution is rather close in en- 
ergy a static stripe phase is quite unlikely to exist in 
Tl- (Hg-) based cuprates. It is rather more meaningful 
to have in mind a picture where stripes appear as snap- 
shots in the dynamical evolution of the system. In this 
context it is interesting to observe that for all dopings 
diagonal stripes are higher in energy than the polaron 
structure. The low doping phase may resemble that of a 
'stripe nematic— where orientational symmetry is bro- 
ken but quantum fluctuations mainly lead to a restora- 
tion of translational symmetry. 

IV. CONCLUSIONS 

Within the inhomogeneous GA we have investigated 
metallic stripe states in the one-band extended Hubbard 
model for a fixed value of the on-site repulsion U/t = 10 
as appropriate for the cuprates. Metallic stripes were 
shown to be favored by a negative t'/t and the physics 
of this behavior has been clarified in terms of the com- 
petition between the kinetic energy along the stripe and 



FIG. 8: Energies per hole for vertical (a) and diagonal (b) 
stripes as a function of doping, t'/t = —0.4 as appropriate 
for Bi2212. The inset displays the incommensurability versus 
doping as deduced from (a). 



the creation energy of a stripe. In this regard the HF ap- 
proximation has been shown to lead to gross quantitative 
and often qualitative errors. 

Generally the present results are in agreement with 
our previous findings in the 3-band model, in particular 
the behavior of the incommensurability and the chemi- 
cal potential, showing that our previous results in the 3- 
band model are quite robust. On the other hand there 
are also important differences. In the present model SC 
stripes are more stable at low doping as compared to BC 
stripes. This is in agreement with DMFT computations 37 
and therefore seems to be a robust feature of the one- 
band model and not to depend on the approximation. 

One should keep in mind that effective models, like 
the one-band Hubbard model, are in principle designed 
for the description of electronic excitations in a diluted, 
uniform phase. Usually it is not guaranteed that they will 
give correct results when one compares the relative sta- 
bility of dense non-uniform phases^ Therefore we con- 
sider the 3-band result for the stripe state symmetry as 
more reliable for the description of real cuprates than the 
one-band calculation. Further support to this conclusion 
comes from the analysis of transport datai 28 i 29 Indeed 
within the 3-band Hamiltonian where BC stripes consti- 
tute the most stable low doping state we have explained 24 
an anomalous vanishing of the transport coefficientsS2i2£ 
as due to the particular band structure of BC stripes 
which turns out to be particle-hole symmetric close to 
the Fermi level. Since the band structure of half-filled SC 
stripes is not particlc-holc symmetric around the Fermi 
level (cf. Fig. |2J) the one-band model cannot account for 
the observed vanishing of transport coefficients! 2 ^* 2 ^ 

When the next-nearest neighbor hopping is fixed to 
t'/t = —0.2 stripes turn out to be half-filled at low dop- 
ing. The behavior of the incommensurability and chem- 
ical potential is in good agreement with experiment in 
LSCO. Our investigation therefore supports the LDA 
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based one band-mapping of Ref. ^ where a similar value 
for the next-nearest neighbor hopping has been reported 
for LSCO materials. 

As doping increases the energy becomes independent 
of the charge topology. SC vertical, BC vertical, diago- 
nal stripes and a polaron lattice all become degenerate. 
We associate this with a transition from a state with ori- 
entational order to an isotropic stated This quantum 
melting of the stripe state for x > 0.2 is also reminiscent 
of the existence of a quantum critical point in the slightly 
overdoped regime as has been proposed by Castellani et 

Since stripes are charged the lateral motion of stripes 
is optically active. Therefore the decrease in the energy 
gap between SC and BC stripes should reflect in the soft- 
ening of an optically active electronic mode. We have 
explicitly shown this behavior in a GA+RPA computa- 
tion in the three-band model where the softening of a 
stripe phasonS can explain the shift of the mid-infrared 
peak in optical conductivity experiments on LSCOi^iM 
We expect a similar behavior in the charge response of 



the one-band model. 

Our present results allow for predictions for other 
cuprate families. For larger values of \t'\ we obtain lower 
values of v. Therefore according to Eq. © we predict a 
steeper curve of e vs. x for small x in Tl and Hg based 
compounds. As \t'\ increases we also find that the rela- 
tive stability of stripes with respect to isolated polarons 
decreases. This loss of stability may indicate that stripes 
have a more dynamical character in these materials. 

Having in mind the differences and similarities between 
the one-band and the three-band description one can still 
obtain very useful qualitative information from the one- 
band description. Since within the present work we have 
set up the relevant parameter set and stable stripe saddle- 
points for the LSCO compounds this represents the first 
step for an investigation of the associated charge and 
spin dynamics. Furtheron the one-band model allows for 
the investigation of larger clusters and thus the struc- 
ture of these excitations in q-space can be investigated 
in more detail than possible within the 3-band Hamilto- 
nian. Work in this direction is in progress. 
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